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Abstract 

Some aspects of direct ion transfer across the water/ 1,2- dichloroethane 
are analyzed using a very simple model based on thermodynamic con- 
siderations. It was concluded that ion solvation by water molecules 
may occur in some particular cases in the organic phase, delivering 
an important contribution to the Gibbs free energy of ion transfer be- 
tween the aqueous and the organic phase. In general terms, this par- 
ticular type of transfer should be favored in the case of highly charged 
small ions at interfaces with a relatively low surface tension and a 
large difference between the reciprocal of the corresponding dielectric 
constants. 
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1 Introduction 



The understanding of charge transfer processes across liquid / liquid inter- 
faces is particularly important in many areas of chemistry, physical chemistry 
and biology. In the last decades a variety of electrochemical methods have 
been adapted to study these processes. Cyclic voltammetry, chronopoten- 
tiometry, impedance spectroscopy applied to macro and more recently to 
microinterfaces between two immiscible liquids are currently being used 0- 
[||]. Nevertheless, interpretation of the experimental data requires in many 
cases a proper knowledge of the interface HQ. In recent years new exper- 
imental techniques have been developed, which allow to obtain molecular 
level information concerning the interface. Second harmonic generation at 
the liquid / liquid interface provides information about the state of molecules, 
especially the surfactants adsorbed at the interface pf] [p[ fLOf . Scanning elec- 



trochemical microscopy studies of charge transfer processes across ITIES are 
in progress The concentration profiles on both sides of the phase 
boundary can be probed and the thickness of the mixed solvent layer can 
be evaluated by scanning the ultramicroelectrode tip|TB[. More recently, IR 
reflection spectroscopy [|14]] and specular reflection of neutrons Jl5]| have been 
applied yielding information at the molecular level. Theoretical calculations 
should provide an important contribution to give insight into interfacial struc- 
ture and the mechanism of transfer across ITIES; several authors have used 
molecular dynamics to model the interface, the ion transfer |T6| [17 and the 
water /oil interface in the presence of micelles fl8|] or in the presence of a mono- 
layer of amphipatic molecules|T^]. Benjamin |^U| has employed equilibrium 
free energy calculations and equilibrium dynamics trajectories to provide a 
detailed microscopic picture of the ion transfer processes. He showed that 
the transfer into the organic phase is accompanied by a hydration shell of 
water molecules and may also involve ion-pairing. A lattice gas model has 
been recently used by Pereira et al. [^T|and Schmickler [22] to clarify a num- 
ber of experimental features of liquid-liquid interfaces, including capacity 
and ion transfer. On the other hand, continuum models have also been em- 
ployed to understand ion transfer between immiscible liquids. A review on 
this topic has been given by Markin and Volkov|23|j, where the Born model, 
non-linear dielectric effects and the use of non-local electrostatic methods 
were analyzed. Although quantitatively poor, the Born theory was found to 
provide a qualitative understanding of ion transfer processes from an aqueous 
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to an organic phase, concerning the effect of static permittivity and ion size. 
Among continuum models, it is also worth mentioning the contribution of 
Abraham et al. in a sequence of papers pi| P5fl []27| P5| , where calculations 
of free energies and entropies of solvation were performed for a number of 
systems. We shall refer to these papers as A&Ln, the number n denoting 
the number with which each part was originally labelled. 

It is the purpose of the present work to revive and shear new light on 
one of the aspects considered in ASzL5, which is that of ion hydration in 
the organic phase. We add a new feature to their model that introduces an 
extra term to the free energy of transfer which was previously ignored . This 
accounts for the interaction between the water layer surrounding the ion in 
the organic phase and the organic solvent. We also use the Born model to 
discuss some qualitative features of the present formulation. 



2 The model 

We consider the direct transfer of alkaline and alkaline earth ions, indicating 
this process as: 

M +Z (aq) # M +z (org) (1) 



As pointed out in A&L5|28fl in their calculation of partition coefficients, 
the free energy of transfer AG 1 of the ion from water to the organic solvent 
contains in principle neutral and electrostatic terms. Since the former con- 
tributes to AG* only to a minor extent, we shall ignore it in most of the 
discussion given below. 

Table I shows the Standard Gibbs energy AG® r of this reaction as given 
by Sabela et al. compared with the corresponding values for several 

anions. In contrast to the transfer of anions, where large differences are ob- 
served, the transfer energy for the smaller alkaline ions appears to be rather 
size-insensitive. An interpretation of these results could be attempted in 
terms of a continuum model of the solvent with inclusion of non linear di- 
electric effects or non local electrostatic models, as proposed in an interesting 



review by Markin and Volkov p3| . However, we first want to draw attention 
here to an alternative model, which is also based on a continuum model for 
the solvent, but introduces a feature that has been absent in most of the 
previous theoretical considerations. This characteristic is the allowance for a 
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remaining solvation sheet belonging to the aqueous phase even after the ion 
has been transferred to the organic phase. In fact, although Abraham and 
Liszi did consider the existence of a hydration sheet in the organic phase, 



a relevant contribution to the free energy of transfer was omitted there, as 
discussed below. At the very preliminary stage of development of these spec- 
ulations, we will not seek a close agreement with the experiment, but rather 
point out under which circumstances this kind of events may occur. More 
sophisticated theoretical considerations and further experimental evidence 
will allow the confirmation or rebuttal of the present ideas. 

We turn now to consider a naive but enlightening model for the trans- 
fer of ions from the aqueous to the organic phase and we later considered 
some improvements to it within the continuum approach. In order to avoid 
confusion, we want to state here what we mean when we write below the 
" radius of the ion" . In the case where we use ideas related to the work of 
Abraham and Liszi, we shall be talking about the crystallographic radius. 
On the other hand, when we apply the concept developed by Rashin and 



Honig , we shall be actually referring to the cavity that the ion produces 
into the solvent. A detailed discussion about this point has been given by 



Rashin and Honig [ 31 1. Our first simplified model is shown in Figure 1 and 
follows the ideas suggested in A&iLb. We assume that the alkaline ion of 
radius r Q eventually may carry with it a solvation sheet of radius r s — r a , so 
that after the transfer it is immersed in a mixed dielectric medium, consti- 
tuted by a sphere of dielectric constant Si, the same as that of the aqueous 
phase, which is in turn embedded in a dielectric continuum of constant 62- 
We shall refer to this kind of ionic transfer as a transfer through a "water 
peel (WP)" mechanism. If we calculate the electrostatic free energy change 
for this process according to the spirit of Born's model, we get: 



2/0/ : 2 1 . 1 



(2) 



where q represents charge of the ion. However, a very important element is 
missing in this model, as can be understood from the following argument. 
The presence of the aqueous solvation sheet makes the ion more stable in 
the organic phase, and this stability increases as r s increases. Thus, on 
electrostatic grounds r s should increase indefinitely to stabilize the ion as 
much as possible. The reason for this not to happen is that the creation of 
the water layer in the organic phase must have some cost, this depending on 
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the radius r s . Note that these considerations are absent in A&L5, where the 
water layer is somewhat arbitrarily (but very reasonably) assumed to be of 
one water molecular diameter. According to the previous argument, we shall 
write the free energy of transfer according to: 



AC ' 

^transf 



2r e 1 



2 l ei 



1 

r 



+ 



+ 4%^r 



(3) 



where 7 denotes the surface tension of a water /organic solvent interface. 
Thus, the three terms on the r.h.s. of eqn. |3] represent the solvation energy 
of the ion in the aqueous phase, the solvation energy in the mixed phase and 
the free energy required to create the new surface originated by the water 
sphere of radius r s , respectively. The last term looks like the solvophobic 



term discussed by Markin and Volkov when considering the contribution 
of the solvophobic effect to the resolvation energy. However, it is clear that 
it has here quite a different meaning, since r s is no longer the ionic radius 
but that of the hydrated ion. 

On the other hand, the transfer of an alkaline ion from the aqueous to 
the organic phase without the aqueous solvation sheet would involve the 
contribution: 



trans 



2"^ 



1 

£2 



1 



(4) 



We are aware that we actually should use instead of this equation a similar 
one taking into account that the cavity radius produced by a particular ion 



will vary in different solvents pi] , but we postpone this discussion for below. 

Thus, the fact whether the ion will be transferred with an aqueous sol- 
vation sheet or not, will depend on the relative magnitudes of the processes 
involved in || and |j. Simplification of eqn. || yields: 



AC ' 

trans 



1 1 

Note that according to this expression AG' tran 



q £ 1 



(5) 



is independent of the 
ionic radius, in agreement with the experimental observation for the smaller 
alkaline ions. We will now analyze under which conditions the processes 
related to |4] or [5] should prevail. We illustrate this point in Figure 2, where 
we make a qualitative plot of AG tra ns and AG' trans as a function of the 
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radius of the aqueous sphere surrounding the ion. Three different situations 
may appear. In the first (Figure 2a) the curves for AG' trans a.nd AG tr(OTS 
intersect in two points, and there exists an infinity of values of r s which 
allow the inequality AG' trans < AG trans . Alkaline ion WP transfer should be 
the predominant mechanism for the ion transfer in this case. In the second 
(Figure 2b), there is only one point in which we get AG' t 



trans 



AG 



trans 



In 

the third case, we get AG' trans > AG trans for all r s , so that ion WP transfer 
should not occur. 

The relevance of the physical parameters involved in the present model 
may be envisaged from the following analysis. Let us set the necessary con- 
dition for a minimum in eqn. |3| according to: 



dAG 



Of 

trans 



q 2 1 


■ 1 


1 


2r 2 


.e 2 





+ 87T7r s 







(6) 



which leads us to the value for the radius of the aqueous sphere at the mini- 
mum: 



1 

c 2 



1 



Substituting r™ m back into eqn. |^, we get: 



(1/3) 



(7) 



a /~iOl min 
^ trans 



3q 2 



1 

So 



1 



(8) 

2 2 rf m - 

J. Thus, the existence of a minimum for 
AG' trans located below AG tra ns (the case depicted shown in Figure 2a) will 
imply the condition : 



which resembles in shape eqn. 



If we now substitute r™ n into eqn. 
inequality: 



(9) 

and rearrange, we arrive at the 



q 2 54-7T 



1 

£2 



1 



This is a very simple mathematical formulation emphasizing the relevance 
of the physical magnitudes which will allow a WP transfer of ions from an 
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aqueous phase to an organic one: the surface tension of the corresponding 
water/organic interface, the ionic radius and the charge of the ion. Thus, for a 
given water/organic interface, the ionic transfer through the WP mechanism 
will be favored for interfaces with the lower surface tension and ions with the 
smaller radii and the higher charges. In order to make a quantitative test of 
these ideas, we have plotted in Figure 3 the standard Gibbs energies for the 
transfer of some univalent ions from water to 1,2-dichloroethane, obtained 



from electrochemical measurements [ 29 j| 30] as a function of the reciprocal of 



the ionic radius. In Figure 3a we used the ionic radius as usually reported in 
textbooks [0], and in Figure 3b we employed the corrected radius as quoted 



by Rashin and Honig |pl| . This corresponds to the choice of ionic radii for 
anions and covalent radii for cations, all increased by 7 %. Although there 
remain some differences between anions and cations, the latter choice seems 
to be more adequate for this representation . In the case of all the anions and 
the bigger cations, an increase of AG tr ans with l/r is found, a fact which was 
predicted by eqn. The differences observed between cations and anions 
may be attributed to the fact of using the same r irrespective of the solvent. 
It has been pointed out in the literature that different radii for the cavity 
surrounding the ions are expected, especially when the solvents involved have 
a different polarity H31| . On the other hand, a remarkable deviation from the 



linear behavior is found in the case of the smaller cations Li + and Na + . 
Furthermore, the experimental AGtrans is not found to change with the ionic 
radius, a feature which was anticipated in the discussion given above by eqn. 
||. Due to the simplicity of the present model, only a qualitative description 
of the direct transfer can be pursued here. However, we can attempt an 
improvement in terms of dielectric saturation effects as proposed in A&L3. 
Thus, eq. || can be replaced by: 



OO 



AC ' =i'f- 

^trans 2 J r 2 
ro 



£2(r £i(r) 



dr + Aivyr 2 (10) 



where £i(r) and £2( r ) are now electric field dependent functions as given in 
A&L3. These can be calculated numerically by solving the equation e = 
n 2 + (coth(!|p) - along with E = , where n denotes internal 
refractivity of the liquid, p the number density, /i the dipole moment and 
kT has the usual meaning. The integral can also be performed numerically 
and numerical minimization of |X0| provides an estimation of r™ m and AGJ r ™°. 
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These are given in Table 2 along with the more approximate values stemming 
from Born's theory. These values point towards a hydration sheet close 
to a molecular diameter, as suggested in A&iLb and the value for AG' trans 
reasonably agrees with the experimental finding for Li + and Na + , where we 
expect WP transfer to take place. Since there is theoretical evidence that the 



surface tension of small drops should decrease with increasing curvature [35 



we expect from eq. |10| that the hydration extent in the organic phase should 
be somewhat larger than what we calculate here. However, the minimum 
described by this equation is relatively flat so that important changes in 
the hydration extent should not appreciably change the values of AG' trans we 
report here. A further comment is also relevant concerning the accuracy of 



eq. 110. We expect the values obtained from it for WP transfer to be much 



more accurate that any estimation of transfer free energy that can be made 
for ion transfer without hydration. In fact, if we assume that the structure of 
the hydration sheet remains essentially unaltered in the organic phase, this 
equation reduces to: 



oo 



AC ' =£-'['- 

^trans 2 J r 2 



£2{r £i(r) 



dr + ATT'jr 2 (11) 



the integral being extended over regions relatively far from the ion. This 
would naturally not be the case for an analogous extension of eq. f| to 
consider dielectric saturation, which would result in: 



2 00 1 

AC 

^^trans 2 J r 2 
ro 



e 2 (r) £i(rl 



dr (12) 



The value of integral is here according to our experience extremely sensi- 
tive to assumptions made concerning the structure of the dielectric media in 
the neighborhood of the ion. 

Returning to our discussion above, it seems that hydration in the organic 
phase should only occur for the smaller ions. However, the model given 
above is only adequate to describe a complete solvation layer, and hydration 
by one or a couple of water molecules may occur. A rough estimation of the 
free energy decrease AG repi that accompanies the replacement of part of the 
organic solvation layer by a single water molecule may be made according 
to the following argument. We divide the free energy change related to this 
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hydration process in two parts. One corresponding to the ion-dipole inter- 
actions AGi-di P and another one corresponding to dipole-dipole interactions 
AG dip _ dip . The change related to the ion-dipole interactions can be written 
as: 



A ^o <? 2 o ^ , Pdce <? 2 ~ 

^^ion-dip — 2 2 ^H 2 < COS [ Oh 2 > H 2 2 pDCE < COS 'OdCE > 

n H 2 r i-H 2 PH 2 n DCE r i-DCE 

(13) 

where n labels the refractive index, /i the dipole moment, p the number 
density and < cos 6 > the average orientation of water (H 2 0) and 1,2- 
dichloroethane (DCE) dipoles respectively The ion-dipole distances denoted 
as Ti_ H2 o and t"i—dce can be taken to be equal to the sum of the ionic radius 
plus the corresponding molecular radius. The average orientation can be es- 
timated using Langevin's equation L(x) = coth(x) — 1/x with x = Ep/kT. 
Thus, the first term on the r.h.s of eq. [13] represents the free energy decrease 
when the ion is solvated by a water dipole, and the second term represents 
the corresponding free energy increase due to desolvation of the correspond- 
ing part of the organic solvent. The values obtained from this equation are 
shown in Table 3 for the different ions considered here. As an estimation for 
the AG dip _ dip term we consider the free energy change corresponding to the 
dissolution of water in 1,2-dichloroethane at the saturation concentration, 
that can in principle be obtained from solubility data of H.%0 in 1, 2 — DCE. 



Since these solubility values reported in the literature [36] vary between 0.1 



M and 10 _3 M, we estimated AG dip _ dip to be in the range 10 — 30KJ/mol. 
Whatever the accurate value, we see that within the present approximation 
practically all the ions considered in Table 3 should prefer solvation with one 
water molecule to be anhydrous in the organic solvent. Moreover, since the 
resulting values for AG® epl = AG° ion _ dip + AG dip _ dip are of the same order 
of magnitude that a reasonable calculation or free energy of transfer with- 
out hydration would yield, we expect that even solvation with a single water 
molecule should deliver a meaningful contribution to this process. 

A more quantitative approach would require a simulation based on realis- 
tic potentials describing the ion- water, ion- 1,2-DCE and all the intermolec- 
ular interactions. In this respect, it is worth mentioning here the interesting 
results of the simulation performed by Benjamin [ 20 1 concerning the transfer 
of an ion from a nonpolar to a polar phase. In the case where the liquid-liquid 



9 



interface is not forced to be sharp but a wider transition region is allowed, an 
important part of the polar solvent solvation shell is already formed around 
the ion when it is still on the nonpolar side of the liquid-liquid interface. It 
would be thus worthwhile to perform this type of simulation with more spe- 
cific potentials for the different ions. This would provide an important test 
to the present ideas. Determination of standard free energies for the transfer 
of alkaline-earth ions will be also helpful to distinguish between a WP or a 
"dry " transfer . In fact, any electrostatic estimation of the standard Gibbs 
energies of the dry transfer should scale in a first approximation with q 2 , as 
predicted by Born's equation. On the other hand, replacing eqn. [?| into eqn. 
H yields: 



AGLr = ^ (4/3) (27vr) (1/3) 



1 i n (2/3) 

(14) 



-£2 £i- 

Thus, the WP transfer free energies are expected to scale with g( 4//3 \ 
which is much weaker than the q 2 dependence given above. 

We can summarize the present results saying that we have clarified some 
aspects of direct ion transfer across the water/1,2 dichloroethane by using a 
very simple model based on thermodynamic considerations. We believe that 
ion solvation by water molecules may occur in some particular cases in the 
organic phase, delivering an important contribution to the Gibbs free energy 
of ion transfer between the aqueous and the organic phase. This particular 
type of transfer should be favored in the case of highly charged small ions, at 
interfaces with a relatively low surface tension and a large difference between 
the reciprocal of the corresponding dielectric constants. 
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Table 1. Standard Gibbs transfer energies of univalent ions from water 
to 1,2-dichloroethane (as quoted in Table 2 of ref. p9|). 



Ion 


AG° tr /kJ mol- 1 


H+ 


53 


Li+ 


57 


Na+ 


57 


K+ 


50 


Rb+ 


42 


Cs+ 


35 


cr 


51 


Br- 


43 


I- 


33 



Table 2. Calculated radii of solvated ion and standard Gibbs energies for 
ion transfer surrounded by a layer of water. The values labelled with (Born) 
were calculated using Born's solvation model while the values labelled with 
(Sat) were calculated taking into account dielectric saturation through eqn. 
|T0| . The parameters employed in the numerical minimization of eq.[l(] were: 
n H2 o = 1.3329, n DCE = 1.4448 , e H2 o(bulk) = 78.54, e DCE (bulk) = 10.23, 
1h 2 o/dce = 25.5 dyn/cm, 



Valence / Z 



1 

2 



(Born) 

2.48 

3.94 



f ^ 

(Sat) 

3.44 

5.22 



AG'ZZ/kJ mot 

(Born) 

35.6 

89.7 



AG' t ^JU mol 

(Sat) 

45.6 

108.8 



Table 3. Ion-dipole contribution to the free energy associated with the 
replacement of part of the organic solvation layer of an ion by a single water 
molecule . All values are in kJ mol^ 1 . 
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Ion 


AG ion - di p/kJ mol' 1 


Li+ 


-117.4 


Na+ 


-91.1 


K+ 


-69.0 


Rb+ 


-62.4 


Cs+ 


-54.5 


F" 


-69.0 


cr 


-49.9 




-45.5 


r 


-39.5 



Figure captions 

Figure 1: Model for the direct transfer of an alkaline ion of radius r$ through 
the water/organic interface. We allow for the possibility that the ion 
carries with it a solvation sheet which is essentially similar to that 
existing in the aqueous phase. The conditions under which this may 
occur are discussed in the text. 



Figure 2: Qualitative plot of the transfer free energy of an ion with an 



aqueous solvation sheet (AG' trans ) , as a function of the radius r s of 
a sphere containing the ion+aqueous layer, a) The horizontal line 
labelled (AG trans ) indicates the transfer free energy of an ion with- 
out the aqueous layer. These schemes correspond to the conditions: 
a) 



«2 



< 



1 

54vr 



£2 



£1 



i 

54tt 



£1 



> 



i 

54?r 



£2 



Figure 3: Electrochemical standard free energies for the transfer of uni- 
valent ions from water to 1 ,2-dichloroethane. Empty circles indicate 
cations, while filled circles denote anions. The reciprocal radii were 
calculated from a)ionic radii as reported in b) corrected radii as 
reported by Rashin and Honig HW. 
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